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bation  of  the  uniformity  of  the  steps. 
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The  two-parameter  class  of  all  40-contractive  two-step  second-order  formulas  was 
derived  for  arbitrary  step  ratios.  Similarly,  the  one-parameter  class  of  all  4 -contractive 
two-step  second-order  formulas  was  derived  for  arbitrary  step  ratios.  A  specific  4  -contractive 
formula  was  derived  which  minimizes  a  measure  of  the  global  truncation  error.  The  one-leg 
implementation  of  the  A  -contractive  formulas  provides  4-stable  methods  for  x  =  \(/)x  for  any 
X(/),  Re  \(/)  <  0  and  any  step  sequence  {An}.  It  was  shown  that,  among  all  4-stable 
two-step  second-order  formulas  for  uniform  steps,  the  4 -contractive  formulas  are  the  only 
ones  for  which  4 -stability  is  conserved  under  any  sufficiently  small  perturbation  of  the 
uniformity  of  the  steps. 

An  efficient  numerical  study,  as  well  as  perturbation  analysis,  of  the  solutions  to  the 
two-junction  interferometer  was  carried  out.  We  obtained  the  dependence  of  the  resonant 
current  on  the  interferometer's  material  parameters  as  well  as  the  previously  undetected 
subharmonic  resonances.  Also  the  shapes  of  the  phase  functions,  both  near  and  far  from 
resonance,  were  found. 
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1.  RESEARCH  OBJECTIVES 


I.  Derive  reliable  and  efficient  variable-step  methods  for  solving  variable  coefficient  and 
nonlinear  differential  equations. 

II.  Use  analytical  and  numerical  methods  for  determining  the  responses  of  interferometers 
and  other  devices  involving  Josephson  junctions. 
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2.  STATUS  OF  RESEARCH  EFFORT 

I.  Stiff  Equations 
a.  Background 

n 

On  a  non-uniform  grid  tn  =*  2  hv%  n  =  1 ,  2,  ...»  hn  =  tn—tn_x%  the  linear  multistep 

**=  i 

formulas  written  in  operator  form  are 

Pnxn~"hnan*n  =  <*) 

A  A 

where  pnxn:=ZaJnx „_k+J,  onxn:=2pn X„_k+J.  For  uniform  steps,  px„  -  E~kp(E)x„ 
and  axn  *  £"/c<7(£)irtf  where  Exn  :=  jcrt+1  and  p(f)  =  2  ayfy,  o(f)  =  S  Pyfy  are  the  familiar 
polynomials.  The  formula  (1)  is  normalized  by 

•.i-  (2) 

7-0 

In  order  to  solve  the  non-linear  system  x  =  /(/,  x),  the  formula  (1)  can  be  implemented  either 
as  a  linear  multistep  (MS)  method. 


Pnxn  xn ) 


(3) 


or  as  a  one-leg  (OL)  method 


Pnxn~"hflf{ontni  OnXn)  ®  0.  (4) 

The  familiar  MS-method  (3)  is  converted  into  its  OL-"twin”  (4)  by  a  permutation  of  /  and  a. 

The  uniform  step  (variable  step)  version  of  (1)  is  said  to  be  stable  with  respect  (w.r.) 
to  Jr  »  Ax  (Jr  »  \(/)x)  at  q  *  h\  ( qn  *  /*„*(/„))  if  the  characteristic  polynomial  p(f)-qo(J) 
(p„(f  )-<?„<*„(£))  satisfies  the  ’’root  condition”.  For  uniform  steps,  this  is  equivalent  to 
boundedness  of  all  solutions  {x„{  of  the  difference  equation;  for  nonuniform  grids  it  is  a 
formal  algebraic  constraint.  The  set  S  of  all  q's  at  which  the  formula  is  stable  is  the  stability 
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region.  The  formula  (l)  is  said  to  be  contractive  *[8]  at  q(qn)  w.r.  to  a  given  norm  ||  •  II  if 
ll*JI  <  II*/i,lL  n  =*  1,  2,  ...»  where  Xn  =  ^n-^+2*  — •  *«)■  The  contractivity 

region  ^  y  #  jj  »s  the  set  of  all  q's  at  which  the  formula  is  contractive  w.r.  to  ||  •  || .  The  formula 
is  said  to  be  /1-stable  (A -contractive  )  if  C ||  #  ^ );  here  C_  denotes  the  closed  left 
half  of  the  //-plane.  The  formula  is  said  to  be  ^40-stable  (/l0-contractive )  if 
R_cS(R  clK  jj  #  || ),  where  R _  denotes  the  closed  negative  real  axis  of  the  //-plane. 

b.  /40-contractivity  results. 

We  derived  the  two-parameter  class  of  all  two-step  ( k  =*  2)  second-orde**  ( p  =  2) 
formulas  which  are  i40-contractive  w.r.  to  the  maximum  norm  II  •  11^  for  arbitrary  step  ratios 
r  «  rn  =  hn/hn_ x  [11].  This  generalizes  the  corresponding  results  for  uniform  steps  given  in 
[8J.  For  the  test  problem  x  *=  A(/)x  with  arbitrary  A(/)  <  0,  the  OL-implementation  of  any 
variable-step  /!0-con tractive  formula  provides  an  40-stable  method  for  arbitrary  step  se¬ 
quences  {hrt}. 


c.  /1-contractivity  results. 


We  also  derived  the  one-leg  parameter  class  of  all  p i  *  k  —  2-^rmuIas  which  are 
^-contractive  w.r.  to  ))  •  )\  x  for  arbitrary  step  ratios  [111.  They  are  defined  by 


«°.n-— O-l-).  A,.,,  - 

«!.»,“  — t *  — r(I  — v)J,  0,  „  = 
1 — 1  “ v)’ 


l+r 
I  —r 
2 
1 

l+r 


+  Kr~2 1 
2(1  +r) 
v+|v2, 

+  ^1- 

2(1 +r) 


V— 


- V, 

2(  1  +  r) 

— i— V2, 
2(1  -Hr) 


(5) 


where  the  parameter  v  satisfies 


Ini 

r 


£  v  £  1. 


(6) 


♦Numbers  in  brackets  refer  to  the  cumulative  list  of  publications.  Section  3  of  this  report. 


For  x  *  X(/).v  with  arbitrary  X(/),  Re  X(/)  <  0,  the  OL-implemcntation  of  any 
variable-step  /1-contractive  formula  provides  an  ^-stable  method  for  arbitrary  step  sequences 
{/ij.  We  proved  [11]  that  any  formula  with  p  =  k  *  2  is  A  -contractive  w.r.  to  II  •  II  *  iff  it  is 
G-stable,  i.e.  /t-contractive  w.r  to  some  G-norm  (a  norm  associated  with  a  positive  definite 
quadratic  vector  form).  For  any  given  /t-contractive  p  =  k  =  2  formula  with  r  =  1  (uniform 
steps),  specified  by  a  parameter  value  v  =  Vj,  0  <  v,  <  1,  ^-contractive  extensions  to 
non-uniform  steps  (r&  1)  can  be  defined  in  such  a  way  as  to  keep  the  G-norm  fixed  w.r.  to  n 
[11].  If  this  is  done  then,  for  an  arbitrary  step  sequence  {An},  G-stability  is  assured  for  any 
dissipative  (monotone  negative)  nonlinear  system. 

d.  Optimal  A -contractive  method. 

An  "optimally  accurate"  /40-contractive  p  =  k  =  2  formula  for  uniform  steps  was 
obtained  in  [9]  by  minimizing  a  bound  of  the  global  OL-error.  By  minimizing  the  same 
objective  function  over  the  class  of  all  A  -contractive  p  =  k  =  2-formulas  with  r  =  1,  we  found 
[8]  the  "optimal"  A  -contractive  formula 


+  Txn-h(l'xn_2  +  |i„_,  +  f*„>  =  0 


(7) 


which  is  associated  with  v  =  v,  *  j.  This  formula  was  generalized  [11]  by  extending  its 
minimality  property  to  the  variable  step  case.  It  is  defined  by  v  =  2r/(2r  +  l)  and  has 
coefficients 


“I 

a2‘ 


(l+r)(l+2r) 
1  +r 

“l+2  r  , 

1  +2r+2r 

(!+/•)(  I  +2r). 


.  Po  - 


r(l+r) 


(l+2r) 

=  Po 


Pi  ’ 


\  +  2r+2T 
(1+2  r)2 


(8) 
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e.  ^-stability  for  variable  steps. 

Given  any  p  m  k  =  2-formula  which  is  /4-stable  for  r  «  1,  we  studied  whether 
^-stability  will  be  preserved  under  any  perturbation  of  r  around  r  —  1?  We  proved  [1]  that 
for  r  =  l  +  €  the  ^-stability  constraints  remain  satisfied  to  first  order  in  |  e  |  iff  the  formula  is 
/4-contractive.  Any  other  /4-stable  formula  can  be  destabilized  for  certain  \(/),  either  with 
increasing  steps  (e  >  0)  (as  is  the  case  for  the  familiar  two-step  second-order  backward 
diffentiation  formula  used  in  many  of  the  popular  software  packages)  or  decreasing  steps 
(e  <  0).  A  realization  of  this  result  using  geometric  step  sequences  hn  =  h{rn~l  and  A(o/n) 
=  Ajr1-”  is  given  in  [11].  In  this  case,  the  OL-difference  Equation  has  constant  coefficients 
and  it  solutions  remain  bounded  for  n  «  iff  the  algebraic  /4-stability  constraints  are 


satisfied. 
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II.  Josephson  Intefcrometer  Equations 


a.  Background 

A  Josephson  junction  consists  of  two  superconductors  separated  by  an  extremely  thin 
dielectric  barrier.  The  order  parameter  <f>  (phase  difference  between  the  wave  functions  of  the 
tow  superconductors)  satisfies  the  well-known  Josephson  relations.  In  one  space  dimension, 
these  relations,  when  combined  with  Maxwell’s  equations,  lead  to  the  equation 


<92<f> 

1? 


d(f>  , 

—  =  k  sin  <f>. 


(9) 


For  the  purpose  of  describing  the  Josephson  ac-effect,  one  is  seeking  running  solutions  of  (9) 

i.e.,  solutions  <*>  such  that  d<>/d/  is  periodic  of  a  certain  period  (2 ir/w).  In  a  "voltage-driven" 

junction,  the  frequency  to  -  which  is  proportional  to  the  voltage  -  is  prescribed  and  the  current 

is  given  by  the  difference  in  at  the  two  boundaries  of  the  junction.  In  a  "current-driven" 

junction,  one  has  to  find  #  as  well  as  its  associated  voltage  (frequency)  <o.  In  point  junctions, 
2 

the  term  — -  in  (9)  is  negligible.  We  have  been  interested  in  studying  properties  of  interfer- 
Ox2 

ometers  which  consist  of  a  number  of  point  junctions  connected  in  various  manners  and  loaded 
by  currents. 

b.  Interferometer  results 

We  studied  the  symmetric,  two-junction  current  driven  interferometer,  whose  equa¬ 
tions  are 


/?0,  +a<*>,  +  sin  <>,  +  A-,(<>|-02)  =  //2  +  /c 
P<t>2  4*  n<£>2  +  ^ 


(10) 


here  A  is  proportional  to  the  inductance,  !c  is  a  control-current  and  *  <o/  +  periodic  of 
period  2ir/wf  where  is  proportional  to  the  voltage.  The  following  results  were  obtained: 
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i)  There  are  subharmonic  resonances  in  the  /—  V  curve,  occurring  around  one-third  of  the 
resonant  frequency  for  small  A,  and  around  one-half  of  resonance  for  A~l.  These  subharmon¬ 
ics  were  previously  undetected,  but  there  is  now  experimental  evidence  of  their  presence. 


ii)  Near  resonance,  the  periodic  parts  of  the  two  phase  functions  are  approximately  sinusoidal 
and  are  of  opposite  signs.  However,  they  have  an  almost  saw-tooth  shape  for  small  voltages. 


iii)  The  correct  dependence  of  the  resonant  current  on  the  two  intrinsic  parameters  A  and 

T  =  was  found.  Previously,  only  the  approximate  dependence  on  T  was  known. 

v2  a 
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